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After a general introduction to the field, we describe some recent results concerning disorder ef- 
fects on both 'random walk models', where the random walk is a dynamical process generated by 
local transition rules, and on 'polymer models', where each random walk trajectory representing 
the configuration of a polymer chain is associated to a global Boltzmann weight. For random walk 
models, we explain, on the specific examples of the Sinai model and of the trap model, how disorder 
induces anomalous diffusion, aging behaviours and Golosov localization, and how these properties 
can be understood via a strong disorder renormalization approach. For polymer models, we dis- 
cuss the critical properties of various derealization transitions involving random polymers. We 
first summarize some recent progresses in the general theory of random critical points : thermody- 
namic observables are not self-averaging at criticality whenever disorder is relevant, and this lack 
of self-averaging is directly related to the probability distribution of pseudo-critical temperatures 
T c (i,L) over the ensemble of samples (i) of size L. We describe the results of this analysis for the 
bidimensional wetting and for the Poland-Scheraga model of DNA denaturation. 

I. INTRODUCTION 

Random walks and diffusion processes have been the subject of constant interest in mathematics and in physics 
during the last century, for two main reasons : on one hand, they play a central role in probability theory, and present 
a large number of very nice mathematical properties; on the other hand, they naturally appear in a great variety 
of situations in physics and in biology. We refer the interested reader to the various reviews that have appeared in 
2005 to celebrate the 100 years of Einstein's theory of Brownian motion, written either by mathematicians |8Cll87|. 
by physicists 0,ESH3j or by 'biology-oriented' physicists |32ll5^| . 

In the following, we will be interested into the effects of quenched disorder on random walks : are the usual 
properties of random walks stable with respect to the presence of some disorder or inhomogeneity ? if not, what are 
the new properties induced by disorder? 

To answer these questions, it is important to distinguish from the very beginning what we will call 'random walk 
models', where a random walk trajectory is generated by local dynamical rules, from what we will call 'polymer 
models' , where each random walk trajectory is associated to a global Boltzmann weight. 

A. Random walk models : trajectories generated by local dynamical rules 

The usual random walk model is defined on a one-dimensional lattice {n £ Z} by the following local dynamical 
rule : the particle starts at n — at time t — 0, and then at each time step, the particle jumps either on the right 
or on the left with probabilities (1/2, 1/2). The position n(t) reached at time t scales with the well-known diffusive 
behaviour n(t) ~ i 1 / 2 . 

Quenched disorder can be introduced in two different ways in the local dynamical rules, either in the relative 
probabilities to jump on the right or the the left, or in the time spent on a site before jumping out of it : these two 
types define the two basic 'random random walk' models, known respectively as the Sinai model and the trap model. 

1. Sinai model : random bias on each site 



The so c alled 'Sinai model' has interested both the mathematicians since the works of Solomon |107| . Kesten et al. 
[82|. Sinai |l04j . and the physicists since the works of Alexander et aZ.@], Derrida-Pomeau |4jj. There ha ve been m any 
developments in the two communities for the last thir ty y ears: we refer the reader to the recent reviews [TOrlfTl^ for 
the mathematical side and to the reviews 0, |§3> [7^ for the various physicists approaches. 

In the Sinai model, each site n of the lattice is characterized by a quenched random variable uj n G]0, 1[ representing 
the probability to jump on the right, whereas (1 — u„) represents the probability to jump on the left. Once these 
random variables (u;„) have been drawn on the full line, the random random walk is generated as follows : the particle 
starts at n — at time t = 0, and then at each time step, the particle being on site n jumps either on the right or on 
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the left with probabilities (w n , 1 — u n ). Note that this process naturally appears in various contexts, for instance in 
the dynamics of a domain wall in the random field Ising chain |55| or in the unzipping of DNA in the presence of an 
external force 88]. 

2. Trap model : random trapping time on each site 

Trap models propose a very simple mechanism for aging [Tflj . and have been thus much studied from this point 
of view, either in the mean field version H Il9l |23| , or in the one-dimensional version in iii in m mil, which 
presents two characteristic time scales for aging, in contrast with the mean field case. 

In this ID version of the trap model, each site n of the lattice is characterized by a quenched random variable r„ 
representing the mean trapping time spent at site n before jumping either to the right or to the left with probabilities 
(1/2, 1/2). Once these random variables (r„) have been drawn on the full line, the continuous random walk is defined 
by the master equation 

dP t (n) _ P t (n) | Pt{n-l) | P t {n + 1) 
dt r„ 2t„_i 2r„ + i 

Again, this process naturally appears in various physical contexts 0, fill Irs! . 

3. Disorder effects : anomalous diffusion, aging and Golosov localization 

The random random walks such as the Sinai walk and the trap model defined above, represent simple dynamical 
models containing quenched disorder, in which many properties such as aging that exist in more complex systems 0, 
|23,|33 can be studied in details. The first effect of disorder is to slow down the dynamics, because disorder effectively 
generates trapping regions from which it is difficult to escape. For instance, the diffusion becomes logarithmic 
x(t) ~ (mi) 2 in the symmetric Sinai model, or can be algebraic x(t) ~ t a with some exponent a smaller than 1/2 in 
the symmetric trap model, if the mean trapping time diverges. Besides this anomalous diffusion length scale x(t) that 
represents the typical distance travelled during time t, there exists a second important length scale y(t) = Xi(i) —X2(t) 
that represents the distance between two independent particles diffusing in the same disordered sample from the same 
initial condition. In the symmetric Sinai model, as discovered by Golosov, the variable y(t) remains finite with 
probability p — 1 in the limit of infinite time : this is the remarkable Golosov localization phenomenon. More 
generally, in the anomalous diffusion phases of the biased Sinai model and of trap models, it turns out out that 
localization still exists, but is only partial : the variable y(t) remains finite with probability < p < 1 in the limit of 
infinite time. 

In Section [n] we explain how these random walks presenting some localization may be studied via a strong disorder 
analysis, which yields results both for one-time observables and for aging properties. 

B. Polymer models : Boltzmann weight on random walk trajectories 

In 'polymer models', random walk trajectories are not considered as the results of a dynamical process, but represent 
conformations of a 'polymer' chain, on which one defines a probability measure. To see the difference, the most well- 
known example is the Self- Avoiding- Walk (SAW) measure in dimension d > 2 defined as follows : among all the 
random walks of N steps, the weight is uniform for all random walks having no self-crossing, and zero otherwise. It 
is clear that this SAW measure cannot be generated by local dynamical rules. 

In the following, we will be interested into phase transitions of directed polymer models. Here one has to distinguish 
two types of models, namely models where the phase transition already exists in the pure case, and models where the 
phase transition only exists in the disorder case. We now introduce one example of each type. 

1. Bidimensional wetting : disorder effects on critical points 

Wetting transitions are in some sense the simplest phase transitions, since they involve linear systems |49| . Let us 
consider a one-dimensional random walk of increments z(a + 1) — z(a) — ±1, starting at 2(0) = 0. The random walk 
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is constrained to remain in the upper half plane z > 0, but gains an adsorption energy e a if z(a) = 0. More precisely, 
the model is defined by the partition function 

(2) 

rw y i<a<Af y 

In the pure case e a = cq, there exists a phase transition between a localized phase at low temperature, characterized 
by an extensive number of contacts at z — 0, and a delocalized phase at low temperature. The critical temperature is 
simply the point where e^ c€ ° = 2, i.e. where the Boltzmann weight of a contact e^ c£ ° exactly compensate the factor 2 
lost at each contact for entropic reason. At T c , the wall z = is exactly reflexive, whereas for T < T c it is attractive 
and for T > T c it is repulsive. For this model, as well as in other types of depinning transitions, one is interested into 
the effect of disorder in the contact energies e a onto the critical properties of the phase transition. 

In Section IIIII we first summarize some recent progresses in the finite size scaling theory of disordered systems : 
thermodynamic observables are not self- averaging at critical points whenever disorder is relevant, and this lack of self- 
averageness at criticality is directly related to the scaling properties of the probability distribution of pseudo-critical 
temperatures T c (i,L) over the ensemble of samples (i) of size L. We then explain how this framework is very useful 
to characterize various derealization transitions involving random polymers, in particular the bidimensional wetting 
transition introduced above, and the Poland-Scheraga model of DNA denaturation. 



2. Directed polymers in random media : transition towards a disorder-dominated phase 
The model of directed polymer in a 1 + d random medium is defined by the following partition function 

Z N (/3)=J2exp 1/3 ( 3 ) 

RW y l<a<N J 

over d— dimensional random walks r(a), where the independent random energies e(a,r) define the random medium. 
This model has attracted a lot of attention because it is directly related to non-equilibrium properties of growth 
models [66j. Within the field of disordered systems, it is also very interesting on its own because it represents a 
'baby-spin-glass' model H, 0, EE Ell presenting a disorder-dominated phase, where the order parameter is an 
'overlap' : two copies of the polymer in the same disordered sam ple have an extensive number of contacts. This 
localization property, that has been now proven by mathematicians |27l l3 l| , is reminiscent of the Golosov localization 
of random walks in random media. From the point of view of anomalous diffusion however, the two types of models 
are very different, since the disorder-dominated phase of directed polymers in random media are super-diffusive x ~ v= 
where £ > 1/2, because the polymer has more chance to find a best global path by making larger excursion, whereas 
random walks in random media are usually sub-diffusive because disorder generates traps for the dynamics. 

In dimensions d < 2, this disorder-dominated phase exists at all temperatures, whereas in dimensions d > 2, there 
exists a phase transition between this disorder-dominated phase at low temperature and a diffusive phase at high 
temperature [36ll75l| . The critical properties of the phase transition in d = 3 have been studied in [42ll84j via standard 
finite-size scaling on disordered averaged observables. However, as explained above for depinning transitions, random 
critical points can be better understood via the knowledge of probability distribution of pseudo-critical temperatures 
T c (i, L) over the ensemble of samples (i) of size L. We are currently studying this distribution, and preliminary results 
show that the critical point is 'unconventional', with different scalings for the shift and the width of this distribution 
& 

II. RANDOM WALKS IN RANDOM MEDIA 
A. Sinai walk and diffusion in a Brownian potential 

1. Continuous and discrete versions of the model 
The continuous version of the Sinai model corresponds to the Langevin equation P"t| 

^L = -U'{x(t))+r,{t) (4) 
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where rj(t) represents the usual thermal noise 

< r){t)rj(t') >= 2T5(t - t') (5) 
and where U(x) is the quenched Brownian potential 

(U(x)-U(y)) 2 = 2a\x-y\ (6) 

More generally, in the whole paper, thermal averages of observables are denoted by < / >, whereas disorder averages 
are denoted by /. 

In the discrete Sinai model on the ID lattice presented in the Introduction, the particle which is on site i has a 
probability u>i of jumping to the right and a probability (1 — w,-) of jumping to the left. The Ui are independent 
random variables in ]0, 1[. The random walk is recurrent only if lnwj = ln(l — u)i), which corresponds to the absence 
of bias of the random potential U(x) © of the continuous version. 

2. Logarithmic anomalous diffusion 

Th e first important property of the Sinai model is the logarithmic scaling of the typical displacement x ~ (lni) 2 
|104| instead of the usual behaviour x ~ \ft of the pure diffusion. This very slow dynamics can be understood by 
the following argument argument |l7j : the time t[x) necessary to reach the point x > will be dominated by the 
Arrhenius factor of e^ Bx associated to the largest barrier B x which should be passed by thermal activation to go from 
the starting point x = to the point x (This approximation by the Arrhenius factor amounts to apply a saddle-point 
method on the exact expression for the first passage time). In a Brownian potential, the typical behaviour of the 
barrier B x ~ ^fx leads to an Arrhenius time t ~ e@^ 7 which indeed corresponds to the scaling x ~ (hit) 2 after 
inversion. 



3. Golosov localization 

The second important property of the Sinai model is the Golosov localization 64] : the distance between two 
independent particles (i.e. with two independent thermal histories) which diffuse in the same disordered sample, 
remains a finite random variable in the limit of infinite time. 

This remarkable phenomenon means that there are sample-dependent areas that concentrate almost all the proba- 
bility weight, and that thermal fluctuations are completely sub-dominant with respect to disorder. 

4- Strong disorder approach 

In the strong disorder approach of the Sinai model the essential idea is to decompose the process Xu,r)(t), 

representing the position of the random walk generated by the thermal noise rj(t) in the random potential Brownian 
U(x), into a sum of two terms 

X{uM (*) = m {u} (t) + y {UM (t) (7) 

• the process m^u^(t) called the "effective dynamics" depends only on the disorder but not on the thermal noise 
: it represents the most probable position of the particle at the moment t. It simply corresponds to the best local 
minimum of the random potential U(x) that the particle has been able to reach at time t. As the escape over a 
potential barrier F requires an Arrhenius time of order tp = r e l3F ', one can study in detail this effective dynamics by 
using a strong disorder real space renormalization which consists in the iterative decimation of the smallest barriers 
remaining in the system. One then associates to time t the renormalized landscape in which only barriers larger than 
the RG scale T = Tint have been kept. The position m^u^(t) ~ T 2 — (Tint) 2 then corresponds at the bottom of the 
renormalized valley at scale r = Tint which contains the initial condition at t — 0. 

• the process y{u,n}{t) represents the thermal fluctuation with respect to effective dynamics. In the limit of infinite 
time, it remains a finite random variable. This very strong result is the Golosov localization phenomenon |(34| : all the 
particles which diffuse in the same sample starting from the same starting point with different thermal noises r) are 
asymptotically concentrated in the same renormalized valley of minimum m^u} (t). More precisely, if one considers the 
first corrections at large time, the probability that a particle is not in the valley corresponding to effective dynamics 
m {u}(t) is °f order l/(lni), in which case the particle is at a distance of order (Int) 2 from mjjit). These events are 
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thus rare (their probability tends towards zero to large time) but they nevertheless dominate certain observables, such 
as the thermal width Ax 2 (t) ~ < y 2 (t) > ~ (lni) 3 which diverges. 

As a consequence of the Golosov localization, the distribution of the rescaled variable X = dj^fy , with respect to 

the thermal noise 77 in a given sample is asymptotically a Dirac delta distribution S(X — M) where M — 7j^gjp is the 
rescaled variable of the effective dynamics. To compute the averaged diffusion front over the samples (or equivalently 
over the initial conditions), it is then enough to study the distribution of M over the samples |52,|j^|. This leads to 
the Kesten law 



P(X) = LT-_I |X 



p \ cosh^> 



4 (-1)" 

7T ^ 2n+ 1 

n=0 



which is an exact result of mathematicians |63j, |65|, |83j • This example explicitly shows how the strong disorder RG 
allows to obtain asymptotic exact results, and gives confidence in the new results of the method concerning finer 
properties. 

Similarly, the rescaled variable for the energy 

_ U(x(0)) - U(x(t)) ^ U(m(0)) - U(m(t)) 

{Tint) ~ (Tint) 1 1 

is entirely given at large time by the effective dynamics. The reduced variable w has the following limit law as t — > 00: 

D(w) = 9{w < 1) (4 - 2w - Ae~ w ) + 6{w > 1) (2e - 4) e~ w (10) 

This law is continuous, like its derivative at w = 1, but the second derivative is discontinuous at w — 1, which 
can seem surprising at first sight. Indeed, for any finite time, the energy distribution is analytic, and it is only in 
the limit of infinite time that the discontinuity appears for the rescaled variable. It is interesting to note that in 
the recent mathematical work over the return time to the origin after time t, another not-analytical asymptotic 
distribution for a rescaled variable also appears. The joint limit distribution of the position X = x fj!, ^ffl and the 



(Tint) 2 

energy w — ^^(^lnt)^*^ may also be computed in Laplace transform 



5. Aging properties 



The two-time diffusion front P(x, t; x' , |0, 0) presents an aging regime in (mi/lnt'). In the rescaled variables 
X = (x/ln 2 t) and X' = (x'/ln 2 t), the diffusion front is again determined by the effective dynamics. The RG 
procedure allows to compute the joint law of the positions {m(t),m(t w )} at two successive times t > t w [U- In 
particular, this two time diffusion front presents a Dirac delta function S(X — X'), which means that the particle can 
be trapped in a valley from which it cannot escape between t' and t. The weight D(t, t w ) of this delta function thus 
represents the probability of having m(t) = m(t w ) 



6. Distribution of the thermal packet 

The asymptotic distribution of the relative position y = x(t) — m(t) with respect to the effective dynamics m(t) 
corresponds to the Boltzmann distribution in an infinite Brownian valley 

/ e -/3E/i(M)) \ 

P{V) ^ \S™dxe-^ + S™dxe-^j {ui jj2} (12) 

where the average is over two Brownian trajectories {U\ 1 U2} forming an infinite valley This formulation is equivalent 
to the Golosov theorem [(54| . The law can be explicitly computed in Laplace transform in terms of Bessel functions 
|92| . and presents in particular the algebraic decay 

p (y) - 4j2 ( 13 ) 
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This can be understood as follows: whereas the Brownian potential U(y) yields a typical decay of order e _,3 v /s J' for the 
Boltzmann factor, there are rare configurations which return close to U ~ at a long distance y, with a probability 
of order l/(y 3 / 2 ). As expected, the correlation of two independent particles in the same sample computed in j§3 



+00 



CM) = lim 2 / dx[P(x,t\x ,0)P(x + l,t\x ,0)] (14) 

t— *oo 

presents the same algebraic decay in I /I 3 / 2 . 



7. Thermal width and rare events 



The algebraic decay for the limit law of the relative position y implies that the second moment < y 2 > diverges 
at infinite time. To obtain its leading behaviour at large time, it is in fact necessary to take into account the following 
rare events |54| : (a) a renormalized valley can have two minima which are almost degenerated in energy; (b) two 
neighbouring barriers can be almost degenerated; (c) a barrier can be near the decimation threshold (r + e). These 
rare events appear with a weak probability of order 1/T, but they give rise to a splitting of the thermal packet into 
two sub-packets, separated by a long distance of order T 2 . As a consequence, these rare events dominate the thermal 
width |5t| 

< x 2 {t) > - < xft) > 2 oc ^(T 2 ) 2 = T(Tlntf (15) 
This behaviour of the thermal width has been measured numerically (2^| . 



8. Discussion 



The strong disorder approach of the Sinai model thus allows to obtain many explicit results, and some of them 
have now been confirmed by mathematicians, in particular the results concerning the weight of the singular part of 
the two time diffusion front [3^| , the statistics of the returns at the origin of effective dynamics [3l| , or properties of 
eigenvalues and eigenvectors of the Fokker-Planck operator [24[ . 

As a final remark, let us briefly explain how the strong disor der approach is related to general theory of slow 
dynamics based on the idea of metastable states 0, |^, |23, Il08j . The idea is to decompose the dynamics into two 
parts. There are on the one hand fast degrees of freedom, which convergence quickly towards a local quasi-equilibrium 
: they correspond to the " metastable states " . On the other hand, there is a slow out-of-equilibrium dynamics which 
corresponds to the evolution of the metastable states. In this language, the strong disorder description of the Sinai 
random walk can be reformulated as follows: 

• the metastable states at time t are the valleys of the renormalized landscape at scale T = Tint: indeed, the 
walkers who were at t = inside this valley have not been able to escape from this valley before time t. 

• In each renormalized valley, there is a quasi-equilibrium described by a Boltzmann distribution inside the valley. 

• the slow dynamics corresponds to the evolution of the renormalized landscape with the scale V = Tint: some 
metastable states disappear and are absorbed by a neighbour. 



B. Biased Sinai walk and related directed trap model 

1. Anomalous diffusion phase 

The introduction of a constant force Fq into the Langevin equation of the Sinai model is very natural. This biased 
model has also interested mathematicians and physicists for a long time, because it presents a series of dynamic phase 
transitions [1711401 182| in terms of the dimcnsionless parameter /1 = FqT/u. In particular, there exists an anomalous 
diffusion phase for < /x < 1, which is characterized by the asymptotic behaviour 

< x{t) > ~ f (16) 



whereas for fj, > 1, the velocity becomes finite: < x(t) > ~ V(/i)t with V([i) = Fq(1 — 1/u). In the anomalous diffusion 
phase, the exact diffusion front is given in terms of Levy stable distributions |l7l l7ll |82( 
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2. Related directed trap model 



It has been proposed for a long time [TJ llS |42| that the biased Sinai model should be asymptotically equivalent 
to a directed trap model defined by the master equation 

dP t (n) _ P t (n) | P t (n-1) 
dt r n T n -i 

in which the r n are independent random variables distributed with the algebraic law 

9(X) * -r^ (18) 

The anomalous diffusion phase < fi < 1 then corresponds to the case where the averaged trapping time is infinite. 

The directed character of this trap model allows to obtain many exact results, since the particle visit sites only 
once in a fixed order, from left to right. In particular, the diffusion front can be expressed in terms of Levy stable 
laws 0| . The thermal width has been exactly computed in p| 



+oo +oo / . s 3 

< An 2 (t) > = J2 ^Pt{n) - £ n p t MP = ( ^ W (19) 

n=0 n=0 I W \ A 1 / 

where /(/i) is some explicit integral 0- The result (|19|l shows that the the thermal packet is spread over a length of 
order i M . On the other hand, the infinite-time limit of the localization parameter for k = 2 computed in |35| 



2tt 1 - e ie (M+i) 



(20) 



shows that Yi is finite in the full phase < fi < 1 and vanishes in the limit fj, = 1. How can this property coexist with 
the result H19|) for the thermal width ? The numerical simulations of [3{| show that for a single sample at fixed t, the 
probability distribution Pt(n) is made out of a few sharp peaks that have a finite weight but that are at a distance 
of order i M . This explains why at the same time, there is a finite probability to find two particles at the same site at 
infinite time, even if the thermal width diverges as t 2fi at large time. 



3. Strong disorder approach 

The strong disorder RG presented in the previous section for the symmetric Sinai model can be extended to the 
biased case, but the obtained results are exact in the limit of infinite time t — » oo only if the bias is very small /i — > 
[Hi) : for instance, the RG yields an exponential diffusion front for the rescaled variable X = 2^ that coincide with 
the exact result involving a Levy distribution ^3,H3 on ly m the limit [i — > 0. The reason why the effective dynamics 
is not exact any more when /1 is finite, is that the distribution of the barriers against the bias converges towards an 
exponential distribution of finite width proportional to This shows that the localization of the full thermal packet 
in a single renormalizcd valley at large time, which is valid in the limit fi — > 0, is not exact any more for finite \i. It 
is thus necessary to generalize the strong disorder RG approach to include the spreading of the thermal packet into 
several renormalized valleys. This can be done as explained in 94], and from the description of the diffusion front in 
each sample, one can compute exact series expansion in fi for all observables. In particular, the explicit computations 
up to order order /j, 2 [94[ of the diffusion front for the rescaled variable X = ^ , of the thermal width 



Um ^±LW ± = M (21n2) + /i 2 [-— + 2 In 2(ln 2 - 2 + 2 7g )] + 0(^ 3 ) (21) 

t^oo 

and of the localization parameter 

Y 2 (n) = l-//(21n2) + M 2 (41n2— ) + (9( M 3 ) (22) 

o 

coincide with the series expansions of the corresponding results, for the diffusion front |l7l 182^ , for the thermal width 
given in Eq (|19fl and for the localization parameter given in Eq (|20|) . These comparisons with exact results obtained 
independently shows that the generalized RG procedure is exact order by order in //. To compute observables at 
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order /j, n , it is thus enough to consider that the diffusion front is spread over (1 + n) traps and to average over the 
samples with the appropriate measure |94| . This approach thus allows to understand how the anomalous diffusion 
phase < /i < 1 presents at the same time a diverging thermal width as i 2/i (|21|1 together with a finite probability 
(|22[l of finding two particles in the same trap at large time. 
Moreover, the strong disorder approach yields a quantitative mapping between the biased Sinai model and the 
directed trap model at large times. All the results for the directed trap model can then be translated for the biased 
Sinai model, one simply has to replace traps by renormalized valleys. In particular, the thermal width has for 
expansion 



< Ax 2 (t) > _ (a 2 /3 3 ) 



2/< 



t 2n 



a 2 (3 4 



^^ + [-^ + 21n2(ln2-2-2 7£ )]4 + 0(i) 



(23) 



In conclusion, the anomalous diffusion phase x ~ with < /j, < 1 of the biased Sinai model is characterized by 
a localization on several renormalized valleys, whose positions and weights can be described sample by sample. The 
generalized strong disorder approach allows to compute all observables via a systematic series expansions in \x. 



C. Symmetric trap model 

The symmetric trap model is defined by the master equation QJ, where the trapping times r„ = e /3En are defined 
in terms of random energies E n distributed with the following exponential distribution 

p(E) = 9(E)^e~^ (24) 

This choice of exponential distribution comes from the exponential tail of the Gumbel distribution which represents 
an important universality class for extreme statistics. The exponential distribution of energies (|24J) translates for the 
trapping time r = e@ E into the algebraic law 

q (r)=e(r>l)-^ (25) 



with exponent 



M = 7f (26) 



At low temperature T <T g , the average trapping time J drrg(r) diverges, and this directly leads to aging effects. 

In the symmetric model, each site can be visited several times, which leads to an essential change in the strong 
disorder approach j^] with respect to the directed trap model : a trap of the renormalized landscape will be charac- 
terized by two important times, namely (i) its trapping time Tj, which represents the typical time of exit towards its 
immediate neighbours (ii) its escape time, which represents the time needed to reach a deeper trap. 

In the limit /i — > 0, the following effective dynamics becomes exact: at time t, the particle starting from the origin at 
t = will be at time t either on the first renormalized trap M+ at distance 1+ on the right or on the first renormalized 
trap M_ at distance Z_ on the left. The weight of the trap M + is simply the probability /_/ (l + +1-) of reaching M + 
before M_ in a flat landscape 

Peff(x,t)^ r ^j-5(x + l.) + T ^ J -6(x-l + ) (27) 

This shows that the dynamics always remains out-of-equilibrium : the weights of the two traps are not given by 
Boltzmann factors, they do not even depend on the energies of the traps, but only on their distances to the origin. 

In the rescaled variable X = x/£(t) where ~ represents the anomalous diffusion length scale, the 

average of the diffusion front H27|l over the samples reads |9fij 

g,(X) = e-W / due- u — — + 0(/x) (28) 
Jo \X\ + u 

The localization parameters, which represents the averages over the samples of the probabilities to find k independent 
particles on the same site, are given by [95| 



Y k {n) = hm £ P*(n,t|0,0) = + 0( M ) (29) 
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This result is in agreement with the numerical simulations of Bertin and Bouchaud |Tlj who have obtained Yq, — > 2/3 
and Y3 — s- 1/2 in the limit \i — » 0. 
The thermal width reads 



c 2 (M) = t lim „ = l + 0(/x) (30) 

and more generally, the others thermal cumulants can be derived from the generating function 

Z^s) = In < er s Wl > = f + °° d\e~ x \ ( ^ coth ^ - + 0( M ) (31) 



2 2 



The two-particle correlation function reads 



+00 +00 



C(l,t) =J2J2 P(n,t\0,0)P(m,t\0,0)6 l!ln _ ml ^ Y^S^ + 



1 



"(a')) 



-C, 



The weight of the 5 peak at the origin correspond as it should to the localization parameter Y2 = 2/3 + O(fi) l|29|l . 
whereas the second term involves the following scaling function 

C fi (\) = e - x ~ + 0(n) (32) 

Finally, one can also obtain results for the two aging correlations 95], namely for the probability H(t + t w , t w ) of no 

jump during the time interval [t w ,t w + t], that presents a sub-aging scaling form in t/tw \ and for the probability 
C(t + t w , t w ) of being at time (t + t w ) in the trap where it was at time t w , that presents an aging scaling in t/t w . 
The strong disorder approach can also be used to obtain results on the non- linear response in the aging regime [97| • 



III. WETTING AND OTHER DEPINNING TRANSITIONS : RANDOM CRITICAL PROPERTIES 



A. General introduction on random critical points 

1. Harris criterion to determine disorder relevance near 2d pure critical points 

The stability of pure critical points with respect to weak bond disorder is governed by the Harris criterion |67| : 
near a second order phase transition in dimension d, the bond disorder is irrelevant if the correlation length exponent 
vp = fpure > 2/d. On the contrary if vp < 2/d, disorder is relevant and drives the system towards a random fixed 
point characterized by new critical exponents. 

A simple argument to understand Harris criterion is the following. The pure system at a temperature T =/= T c 
is characterized by a correlation length £(T) ~ t~ Vp , where t = \T C — T\ represents the distance to criticality, and 
up the correlation length exponent. The pure system can be divided into nearly independent subsamples of volume 
V ~ £ d (T) ~ t~ dvp . In the presence of an additional weak bond disorder, the averaged bond value (1/V) X^ev h 
seen in a volume V will present fluctuations of order 1 / y/V. So the fluctuations of critical temperatures among the 
volumes of size V will be of order AT C (V) ~ l/VV ~ £,~ d / 2 (T). Disorder will be irrelevant if these fluctuations 
AT C (V) ~ fAvp/i i-, ecomes negligeable with respect to t = \T C — T\ in the limit t — » where the critical point is 
approached. 



2. General bound ufs > 2/d for random systems, and the possible existence of two distinct exponents v 



There exists a general bound for the finite-size correlation length exponent vfs > 2/d in disordered systems 




which essentially means that a random critical point should itself be stable with respect to the addition of disorder, as 
in the Harris criterion argument given above. However, this general bound has to be understood with the subtleties 
explained in |30j ]. In so-called 'conventional' random critical points, there is a single correlation length exponent 
v = i>ps and this single exponent is expected to satisfy the bound. However, there are also 'unconventional' random 
critical points, where there are two different correlation length exponents ! In this case, the 'intrinsic ' correlation 
exponent ^intrinsic can be less than 2/d, whereas the bound holds for the finite-size exponent vps > 2/d. The most 
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well-understood example of the existence of two different correlation length exponents is the random transverse field 
fsing chain (this quantum ID model is equivalent to the 2D classical Ising model with columnar disorder introduced 
by McCoy and Wu (§3), 

which has been studied in great details by D. Fisher via a strong disorder renormalization 
approach [5^ : the exponent v — 1 governs the decay of the typical correlation lnC(r) ~ — r/£, whereas v = 2 governs 
the decay of the averaged correlation ln(C(r)) ~ —?*/£. Exactly at criticality, the typical and averaged correlations 
are also very different, since the typical correlation decays as Ct yp (r) ~ e~ Ws ^, where w is a random variable of order 

1, whereas the averaged correlation is dominated by rare events and decays al gebr aic ally C(r ) ~ l/y ( 3 ~ v ^/ 2 . Other 
critical points with two different correlation length exponents are discussed in [l6|, l5ll llOll Il02l Il05j| . 



3. Lack of self- averaging at random critical points 



In disordered systems, the densities of extensive thermodynamic observables are self-averaging off-criticality, because 
the finiteness of the correlation length £ allows to divide a large sample into independent large sub-samples. At 
criticality however, this 'subdivision' argument breaks down because of th e divergen ce of £ at T c , and a lack of 
self-averaging has been found at criticality whenever disorder is relevant Q, IllOt llll| . More precisely, for a given 
observable X, it is convenient to define its normalized width as 



R X{ T,L)= X ^ L )-^ L )Y (33 ) 
{Xi{T,L)f 



In terms of the correlation length £(T), the following behaviour of Rx(T,L) is expected |ll llllj 

(i) off criticality, when L 3> the system can be divided into nearly independent sub-samples and this leads to 
'Strong Self- Averaging' 

Rx (T, L)~l- off criticality for L > £(T) (34) 

(ii) in the critical region, when L <C £,(T), the system cannot be divided anymore into nearly independent sub- 
samples. In particular at T c (oo) where £ = oo, one can have either 'Weak Self- Averaging' 

Rx{T c (oo),L) ~ — — — for irrelevant disorder (up > 2/d) (35) 

or 'No Self- Averaging' 

Rx(T c (oo), L) ~ Cst for random critical points (36) 

To understand the origin of this lack of self-averaging, it is useful to introduce the notion of sample-dependent 
pseudo-critical temperatures, as we now explain. 



4- Distribution of pseudo-critical temperatures 

Important progresses have been made recently in the understanding of finite size properties of random critical 
points [ll Ya. Hoi I102L flloL lllfj • To each disordered sample (i) of size L, one should first associate a pseudo-critical 
temperature T c (i,L) pi Ho i ITTrl ITTl) . The disorder averaged pseudo-critical critical temperature T® V (L) = T c (i,L) 
satisfies 

T? V (L) - T c (oo) - (37) 

where vp is the correlation length exponent. Eq. I|37() generalizes the analogous relation for pure systems 

TP ure (L) - T c (oo) - L-V up (38) 

The nature of the disordered critical point then depends on the width AT C (L) of the distribution of the pseudo-critical 
temperatures T c (i, L) 



AT C (L) = y/Var[T c (i,L)] - J^(i,L) - (t c (i,L)) (39) 
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When the disorder is irrelevant, the fluctuations of these pseudo-critical temperatures obey the scaling of a central 
limit theorem as in the Harris argument : 

AT C (L) ~ L~ d/2 for irrelevant disorder (40) 

This behaviour was first believed to hold in general 

nouns, 

but was later shown to be wrong in the case of random 
fixed points. In this case, it was argued [U llllf that eq. 1|40|) should be replaced by 

AT C (L) ~ L~ 1Ivr for random critical points (41) 

i.e. the scaling is the same as the L-dependent shift of the averaged pseudo-critical temperature (Eq. ET7|) . The fact 
that these two temperature scales remain the same is then an essential property of random fixed points that leads to 
the lack of self-averaging at criticality. 

Up to now, to our knowledge, the distribution of T r (i, L ) or of another sample-dependent critical parameter has 
been studied for various disordered spin models 0, 0, Il02l , for elastic lines in random media > f° r Poland- 
Scheraga models [9^ and for the selective interface model |l00j . We will describe below the case of Poland-Scheraga 
model, since the variation of one parameter in the Poland-Scheraga model allows to study disorder effects on both 
first order and second order transitions, with either marginal or relevant disorder. 

5. Finite-size scaling in disordered systems 

In pure systems, the finite-size scaling theory relates the critical exponents of the phase transition in the thermo- 
dynamic limit to finite-size effects that can be measured in numerical simulations 25] . In short, this theory says 
that the only important variable is the ratio between the size L of the finite system and the correlation length that 
diverges at the critical point £(T) ~ |T — T c \~ v . So the data Xl(T) for various sizes L should be analysed in terms 
of the appropriate rescaled variable r = (T — T^)L X I V to obtain a master curve of the form L v Xl(T) = 4>(t), where 
y is the exponent describing the decay the observable X exactly at criticality Xl{T c ) ~ 1/L y . Note that using the 
scale-dependent T C (L) instead of the thermodynamic T c = T c (oo) is completely equivalent, since it corresponds to a 
simple translation f = (T — T C {L))L V I V = r + a, as a consequence of (JSSJ. 

In random systems, one has instead data (T) measured at temperature T for various disordered samples (i) of 
size L, and the question is : what is the best way to analyse these data? The usual procedure consists in averaging 
over the samples (i) at fixed (T, L) to apply the pure procedure to these disorder averaged quantities : one tries to 

find a master curve L y X^\T) — <f>(r) in terms of the variable r = (T — T^L 1 ?" . However, this procedure leads to 
extremely large sample-to-sample fluctuations in the critical region, as a consequence of the width of the distribution 
of pseudo-critical temperatures : at a given temperature T, the samples having their pseudo-critical temperature 
T c (i,L) > T are effectively in the low temperature phase, whereas the samples having T c (i,L) < T are effectively 
in the high temperature phase. This mixing of samples in the critical regions makes it very difficult to obtain clean 
results on critical exponents. 

To avoid these difficulties, the following alternative procedure has been proposed [ToL I102L llllj : the data for a 
given sample (i) should be analysed in terms of the rescaling variable r new = (T — T c (i, L))Z 1 ' 1/ with respect to its 
own pseudo-critical temperature T c (i, L). Since T c (i, L) — T® V (L) + ttjAT c (L), where Ui is a random variable of order 
one, and where the mean T£ V {L) and the variance AT C (L) follow the respective behaviours (|37|) and 141|) . it is clear 
that the two procedures are not equivalent, since r new = r + a + is not a simple translation of the constant a as in 
the pure case, because the random variable Ui remains present for arbitrary L. This new way of analysing the data 
allows to reduce very significantly the sample-to-sample fluctuations, as shown for spin models [ToL lOalllll , for the 
non-equilibrium depinning transition of elastic lines in random media , and for disordered polymer models [9jJ . 

B. Application to Poland-Scheraga model with various loop exponent c 

1. Poland-Scheraga model of DNA denaturation 

The Poland-Scheraga model of DNA denaturation |l03l| describes the configuration of the two complementary 
chains as an alternance of bound segments and open loops. Each loop of length I has a polymeric entropic weight 
Af(l) ~ l<t- l /l c , whereas each contact at position a has for Boltzmann weight e _,3eQ . The partition function Z(a;l) 
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with bound-bound ends at (a, 1) thus satisfies the simple recursion relation 

a-i 

Z{a; 1) = e- 0€ * ^ M(a - a')Z(a'; 1) (42) 

a' = l 

The wetting model presented in the Introduction J2J corresponds to a Poland-Scheraga model with loop exponent 
c = 3/2 (this exponent comes from the first return distribution of a ID random walk). 

For DNA denaturation, the question is what is the appropriate value of the loop exponent c? Gaussian loops in 
d = 3 dimensions are characterized by c = d/2 = 3/2. The role of self avoidance within a loop was taken into account 
by Fisher |48(, and yields the bigger value c = dvsAW ~ 1.76, where vsaw is the SAW radius of gyration exponent 
in d = 3. More recently, Kafri et al. [78LF^ | pointed out that the inclusion of the self avoidance of the loop with the 
rest of the chain further increased c to a value c > 2. The transition in the pure model should then be discontinuous, 
as previously found by j2^ in Monte Carlo simulations of SAW 's. The value c ~ 2.11 was in turn measured in three 
dimensional Monte Carlo simulations 0, 0, H(| . 



2. Pure critical properties and disorder relevance : role of exponent c 



In the pure case e a = eo, the model is of course exactly solvable, and it turns out that the critical properties are 
determined by the value of the loop exponent c : for c > 2, the transition is first order with exponent vp = 1, whereas 
for 1 < c < 2 the transition is second order with exponent Vp — l/(c — 1). The Harris criterion thus yields that 
disorder is irrelevant for c < 3/2, marginal for c = 3/2, and relevant for c > 3/2. 

As explained above, the case c = 2.15 where the pure transition is first order is of interest for DNA denaturation. 
The effect of disorder on this transition has been recently debated 

The marginal case c = § has been debated for a long time 0, 1381 I43LI571 IfSflL IsTi llOfll ] and is of special interest since 
it corresponds to two-dimensional wetting as already mentioned. On the analytical side, efforts have focused on the 
small disorder limit : Ref |57| finds a marginally irrelevant disorder where the quenched critical properties are the same 
as in the pure case, up to subdominant logarithmic corrections. Other studies have concluded that that the disorder 
is marginally relevant Isij . On the numerical side, the same debate on the disorder relevance took place. The 

numerical studies of Ref. [HJ and Ref. [38| have concluded that the critical behaviour was indistinguishable from the 
pure transition. On the ot her h and, the numerical study of pointed towards a negative specific heat exponent 
(a < 0), and finally Ref. 109] has been interpreted as an essential singularity in the specific heat, that formally 
corresponds to an exponent a = — oo. 



3. Distribution of pseudo- critical temperatures in Poland-Scheraga models 

In |99j | , the distribution of pseudo-critical temperatures disordered Poland-Scheraga models with different loop expo- 
nents c, corresponding to either a pure second order transition with respectively marginal/relevant disorder according 
to the Harris criterion, or to a pure first-order transition. Two different definitions for the pseudo-critical temperature 
have been considered, one based on the free-energy, and the other based on a sample-replication procedure. Both 
definitions actually yield similar results for critical exponents and scaling distributions, even though they give different 
values for a given sample. This shows that the conclusions that can be obtained from such studies do not depend on 
the precise definition of the pseudo-critical temperature. 

In all cases (c = 1.5, 1.75 and 2.15), and for both definitions of the pseudo-critical temperature the distributions of 
pseudo-critical temperatures were found to follow the scaling form 

Pi(Tc( ^ }) ~ A7UZ) 9 [ X = A7^(Z) ) (43) 

where the scaling distribution g(x) is simply Gaussian 

g(x) = -±=e-* 2 ' 2 (44) 



2tt 

Note however that this Gaus sian distribution is not generic but seems specific to these Poland-Scheraga models, since 
in the selective interface, 1QQJ, the corresponding scaling distribution was found to be very asymmetric. The rescaling 
(I43|l means that the important scalings of the pseudo-critical temperatures distribution are the behaviours of its 
average T® V (L) and width AT C (L) as L varies. 
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For the marginal case c = 3/2 corresponding to two-dimensional wetting, both the width AT C (L) and the shift 
[T c (oo) — T£ V (L)] are found to decay as L~ Uraniom , where the exponent is very close to the pure exponent (v r andom ~ 
2 = Vpure) but disorder is nevertheless relevant since it leads to non self- averaging at criticality. For relevant disorder 
c = 1.75, the width AT C (L) and the shift [T c (oo) — Tf? v (L)] decay with the same new exponent i -1 /"™*™ (where 
Vrandom ~ 2.7 > 2 > ^ pure ) and there is again no self-averaging at criticality. Finally for the value c = 2.15, of 
interest in the context of DNA denaturation, the transition is first-order in the pure case. In the presence of disorder, 
the width AT C (L) ~ L^ 1 / 2 dominates over the shift [T c (oo) — T£ V {L)\ ~ L 1 , i.e. there are two correlation length 
exponents v = 2 and v = 1 that govern respectively the averaged/typical loop distribution. This is reminiscent of 
what happens at strong disorder fixed points [52l u\ , and perhaps indicates that some appropriate strong disorder 
RG description of this random critical point should be possible. 



IV. CONCLUSION 



We have presented some recent results concerning the effects of quenched disorder on random walks and polymer 
models. For random walk models, we have explained how the Golosov localization phenomenon allows to construct a 
strong disorder approach to obtain detailed explicit results on the anomalous diffusion behaviours and on aging prop- 
erties. For disordered polymer models presenting depinning transitions, we have explained how the recent progresses 
in the general theory of random critical points are very useful to understand the critical properties. 

To finish, let us mention some closely related works. In the field of random walks, we have only described here the 
simplest models involving a single random walk in a disordered medium, but models involving many interacting random 
walks in the presence of quenched disorder can be also studied via the strong disorder approach |74| . in particular 
reaction diffusion models HHHK, asymmetric exclusion processes [7^|, contact processes |?3 and zero-range processes 
[77]] • In the field of phase transitions of disordered polymer models, we have restricted our presentation here to the case 
of bidimensional wetting and of Poland-Schcraga model of DNA denaturation, but the idea to study the distribution 
of pseudo-critical temper ature s allows to clarify the critical properties of other equilibrium transitions, such as the 
selective interface model |lOdj | and the directed polymer transition in d = 3 [6lJ mentioned in the introduction |J3}. 
This framework is also very useful for the non-equilibrium depinning transition of elastic lines in random media, where 
the notion of pseudo-critical temperature is replaced by the notion of sample-dependent pseudo-critical force [la . |45| . 
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